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Research Project 


The Interchange No. NCC2-5149 deals with the emerging technology of photonic (or optoelec- 
tronic) integrated circuits (PICs or OEICs). In PICs, optical and electronic components are grown 
together on the same chip. To build such devices and subsystems, one needs to model the entire 
chip. PICs are useful for building components for integrated optical transmitters, integrated optical 
receivers, optical data storage systems, optical interconnects, and optical computers. For example, 
the current commercial rate for optical data transmission is 2.5 gigabits per second, whereas the 
use of shorter pulses to improve optical transmission rates would yield an increase of 400 to 1000 
times. The improved optical data transmitters would be used in telecommunications networks and 
computer local-area networks. Also, these components can be applied to activities in space, such as 
satellite to satellite communications, when the data transmissions are made at optical frequencies. 

The research project consisted of developing accurate computer modeling of electromagnetic 
wave propagation in semiconductors. Such modeling is necessary for the successful development 
of PICs. More specifically, these computer codes would enable the modeling of such devices, 
including their subsystems, such as semiconductor lasers and semiconductor amplifiers in which 
there is femtosecond pulse propagation. Presently, there are no computer codes that could provide 
this modeling. Current codes do not solve the full vector, nonlinear. Maxwell’s equations, which 
are required for these short pulses and also current codes do not solve the semiconductor Bloch 
equations, which are required to accurately describe the material’s interaction with femtosecond 
pulses. The research performed under NCC2-5149 solves the combined Maxwell’s and Bloch’s 
equations. 
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Research Accomplishments 


The main objective of the Joint-Research Interchange NCC2-5149 was to develop computer codes 
for accurate simulation of femtosecond pulse propagation in semiconductor lasers and semicon- 
ductor amplifiers [ 1 ] . The code should take into account all relevant processes such as the interband 
and intraband carrier relaxation mechanisms and the many-body effects arising from the Coulomb 
interaction among charge carriers [2]. This objective was fully accomplished. We made use of a 
previously developed algorithm developed at NASA Ames [3]-[5]. The new algorithm was tested 
on several problems of practical importance. One such problem was related to the amplification 
of femtosecond optical pulses in semiconductors. These results were presented in several interna- 
tional conferences over a period of three years. 

With the help of a postdoctoral fellow, we also investigated the origin of instabilities that can 
lead to the formation of femtosecond pulses in different kinds of lasers. We analyzed the oc- 
currence of absolute instabilities in lasers that contain a dispersive host material with third-order 
nonlinearities. Starting from the Maxwell-Bloch equations, we derived general multimode equa- 
tions to distinguish between convective and absolute instabilities. We find that both self-phase 
modulation and intensity-dependent absorption can dramatically affect the absolute stability of 
such lasers. In particular, the self-pulsing threshold (the so-called second laser threshold) can oc- 
cur at few times the first laser threshold even in good-cavity lasers for which no self-pulsing occurs 
in the absence of intensity-dependent absorption. These results were presented in an international 
conference and published in the form of two papers. 
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An algorithm has been developed that solves the semiconductor Maxwell-Bloch equations [1], without 
making the standard approximations of a slowly-varying envelope (SVEA) and a rotating-wave (RWA). This 
more exact formulation is applied to simulations of the propagation of ultrashort pulses for which the standard 
approximations reach their limits. This development was motivated by the generation of optical pulses as 
short as 8 fs, which has become possible due to recent progress in ultrafast technology. 

Previously, an algorithm was developed for the Maxwell equations [2,3], without making the SVEA 
for calculations of pulse propagation in nonlinear glasses, which exhibit Kerr-like instantaneous nonlinearities. 
It was found that significant differences can occur when the SVEA is not made in Maxwell equations. 
Specifically, light bullets, of 25 fs duration, were found to be stable [4] with the full Maxwell equations, 
whereas previously, calculations with the nonlinear Schroedinger equation had shown them to be unstable 
[5]. More recently [6], using the algorithm developed in references 2 and 3, calculations showed the formation 
of shock waves on the optical carrier wave. Such results are impossible with the SVEA since the carrier wave 
is eliminated from the calculations. 

In this paper, this new algorithm is applied to studies of ultrafast optical pulse propagation in 
nonlinear semiconductor materials, in which many-body effects due to Coulomb interactions are included. In 
addition to the algorithm for Maxwell’s equations, a new algorithm has been developed for the semiconductor 
Bloch equations that does not make the RWA and the two algorithms have been combined into one for the 
coupled semiconductor Maxwell-Bloch equations. In the Bloch equations the relaxation-time approximation 
[1] has been made for the various intraband scattering processes. 

The Maxwell-Bloch equations for pulse propagation in one spatial dimension are the following. 
Assume that the electric field of a pulse that is propagating along the z direction is polarized along the x 
axis and ignore the transverse effects, then the Maxwell equations become 


dD X _ 3Hy 
dt dz 


dH y dE x 

^~dT = “5T 


D x — tO^rEx + Ex i 


where P x is the induced polarization. 

In the case of semiconductors, P x is calculated by using the semiconductor Bloch equations [1], (a 
two band model, one conduction band and one valence band). 
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where n\ and are the occupation probabilities for electron and holes of the wave vector k in the conduction 
and valence bands respectively, Wk — (n£-f-njj -1) is the population inversion, p\ t k and p 2 t k are the dispersive 
and absorptive components of the dipole moment p*, of the wave vector fc, p* — pi } k + *P2,k and the r 
parameters govern various decay processes. The transition energy hujk is varied over a sufficiently large range 
to accurately describe the interaction of an ultrashort optical pulse with the semiconductor. n£(t) and «£(*) 



are determined by first computing the chemical potentials fi e (t) and p h {t) from n*(J) and n£(f) respectively 
and then using the formula for a Fermi-Dirac distribution to find the quasi-equilibrium Fermi distributions 
n£(<) and n£(*) 

The generalized Rabi frequency 17* = fii ( * + *f72,*» the effective transition energy fiA*, which includes 
the band-gap renormalization, and the induced polarization P T are given respectively by 

= + &A t = fc W *-£l^_,|(n; + n£), P x (t) = ^ f pi, k k 2 dk 

n q* k q*k J ° 

Initially an algorithm was developed for the simpler optical Maxwell-Block equations for two-level 
atomic systems [7]. A calculation of self-induced transparency was made for a a 10 fs pulse [8]. Figure 1 
shows the electric field of the pulse at several moments during its propagation inside the medium. The top 
curve on the right side in figure 2 shows the corresponding population inversion at some location as the pulse 
goes by. The top curve on the left side in figure 2 shows the resulting population inversion for self-induced 
transparency when the SVEA and the RWA are made. The remaining curves are a comparison of the two 
methods when the atomic transition frequency is detuned away from the optical carrier frequency. The exact 
method is able to capture the off resonance details that the approximate method is incapable of modeling. 

Next the semiconductor Maxwell-Block equations were solved [8] under the simplifications that the 
Coulomb interaction terms were neglected (the free carrier assumption) and that there were no relaxation 
terms in the equations for the evolution of n e k and n£. Figures 3-6 show colliding pulses, including constructive 
and destructive interference. Figure 7 shows gain curves that were obtained under the free carrier assumption. 

Finally figure 8 shows exciton results. Notice the Is and 2s absorption peaks in the case of the 
2 ps dipole decay time. Here the population inversion was specified at minus one, the generalized Rabi 
frequency was used but the transition energy was not renormalized. The presentation shall include additional 
calculations of propagating and colliding pulses in which all the Coulomb and relaxation terms are included. 

We would like to thank Rolf Binder, Optical Sciences Center, University of Arizona, for his many 
helpful comments concerning algorithm development for the semiconductor Bloch equations. 
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Fig. 1. Two-level atom. Self induced transparency. Fig. 2. Two-level atom. Population inversion. 



Fig. 3. Approaching Pulses 


Fig. 4. Constructive Interference 




Fig. 5. Destructive Interference 


Fig. 6. Separating Pulses 
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We investigate analytically the occurrence of modulation instability in doped fiber lasers and amplifiers using 
a Maxwell-Bloch description for the dopants 1 and without making the usual parabolic-gain approximation. 
We find a new modulation instability occurring near the Rabi frequency, which is not predicted by the conven- 
tional complex Ginzburg-Landau model. We discuss the implications of this new instability for fiber ampli- 
fiers and lasers and analyze the effects of the saturable host absorption on the laser instabilities. Atomic 
detuning is shown to significantly enhance the new modulation instability, in both the normal- and the 
anomalous -dispersion regimes. © 1997 Optical Society of America [S0740-3224(97 >04610-9] 


1. INTRODUCTION 

The onset of instabilities in various kinds of lasers is gen- 
erally studied by use of a rate-equation model based on 
the Maxwell-Bloch equations suitable for a two-level 
atomic system. 1 Such a model, often referred to as the 
(detuned) Lorenz-Haken model, includes atomic polariza- 
tion dynamics and has been used extensively over the 
past two decades. Its use leads to the concept of the sec- 
ond laser threshold, defined as the pump level at which 
the continuous-wave (cw) operation of the laser becomes 
unstable through a Hopf bifurcation, resulting in a self- 
pulsing output. At higher pump levels the laser can en- 
ter into a chaotic regime through a period doubling or an- 
other route to chaos. 1,2 

The advent of fiber lasers during the late 1980’s forces 
several changes to this standard model of laser instabili- 
ties, mainly because the optical fiber, acting as a host to 
the dopants, introduces group- velocity dispersion (GVD) 
and self-phase modulation (SPM), both of which must be 
incorporated for a proper description of the onset of insta- 
bilities in fiber lasers. 3 In fact, these two phenomena 
lead to an instability, known as the modulation instability 
(MI), even in an undoped and unpumped optical fiber. 4 
It is therefore reasonable to expect that the presence of 
GVD and SPM in the host fiber would change the nature 
of instabilities in fiber lasers in comparison with other 
kinds of lasers (gas and solid-state lasers) that are well 
described by the standard Lorenz-Haken model. Apart 
from this propagation-based instability, several other ex- 
planations for the observed instabilities in rare-earth- 
doped fiber lasers have been reported. It was shown 
theoretically and experimentally that the existence of ion 
clusters in heavily Er-doped fiber lasers leads to single- 
mode cw or self-pulsing behavior, whereas the same 
model is also applicable to dual-wavelength or bipolarized 
lasers. 5,6 Other theories and experiments on Er-doped fi- 
ber lasers have shown self-pulsing, chaos, and antiphase 
dynamics between the different polarization eigenstates 
of the optical field. 7,8 The explanation for the self- 
pulsing behavior of Nd-doped fiber lasers has been re- 


ported to be driven by the dynamics of the two field- 
polarization eigenstates that depend on the birefringence 
of the fiber. 9,10 In this paper, however, we focus on the 
propagation-driven MI phenomenon, and do not consider 
any field-polarization dynamics. 

In recent years, the MI phenomenon has been investi- 
gated in doped (active) fibers used to make lasers and 
amplifiers. 3,4,11,12 When doing so, one has to consider the 
nonlinear interaction of the dopants with the optical field. 
A natural choice is to model the dopants as a two-level 
system with an atomic polarization dephasing time T 2 
and a population relaxation time T ] . By far, the most 
popular model employs the parabolic-gain approximation, 
leading to a complex Ginzburg-Landau (CGL) equation 
for the optical field. 3 One study showed that in erbium- 
doped fiber amplifiers, 11 the threshold for MI is consider- 
ably lowered compared with that for undoped fibers. Re- 
cently, Chen et ai 12 included gain dynamics (governed by 
T j) as well as a fast saturable absorber in the model and 
discussed the implications of MI for passively mode- 
locked figure-eight lasers. The full atomic polarization 
dynamics (governed by T 2 ) has, however, been neglected 
so far. 

In this paper we investigate the occurrence of MI be- 
yond the Ginzburg-Landau approximation by considering 
the full T 2 dynamics. We introduce the theoretical 
framework in Section 2 and discuss the consequences for 
amplifiers in Section 3. There, we calculate the steady- 
state solutions and derive a dispersion relation for MI. 
The effect of the population relaxation damping time T x 
and the dipole dephasing time T 2 on MI are studied for 
amplifiers. We find that by cooling the fiber amplifier, 
the bandwidth and the strength of the MI can be greatly 
reduced. In Section 4 we deal with fiber lasers and de- 
rive the dispersion relation for MI at resonance. The role 
of saturable absorption in fiber lasers is investigated, and 
we focus on the possibility of MI occurring in the normal- 
dispersion regime. We find indeed such an instability, 
having its origin in the atomic coherence effects related to 
the atomic polarization dynamics. This new instability is 
found to occur at rather low frequencies (—50 MHz) and 
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may explain the self-starting behavior of mode-locked Nd- 
doped lasers. In Section 5 we discuss the effects of de- 
tuning on the occurrence of MI and discuss the differences 
between normal and anomalous operating regimes. 

2. THEORETICAL FRAMEWORK 

Our starting point is a set of Maxwell -Bloch equations 
that describe the propagation of optical fields in a nonlin- 
ear, dispersive medium doped with two-level atoms (or 
ions). We write the electrical field #(x, y, z, t) and the 
induced material polarization ^(x, y, z, t) as 

%(x, y, 2, t) = 5 xF(x, y)A(z, t) 

x exp[£(/3oZ - w 0 O] + c.c., (1) 

J>(x, y, z , <) = \ xF{x, y)B(z, t) 

x exp[ i ( Pqz - w 0 *)] + c.c., ( 2 ) 

where x is the polarization unit vector of the light as- 
sumed to be linearly polarized along the x axis, Fix, y ) is 
the fiber-mode profile, and (S 0 is the wave number corre- 
sponding to the carrier frequency o> 0 . We assume that 
the field-polarization direction is preserved upon propaga- 
tion and that we are dealing with a polarization- 
preserving single-mode fiber. However, most of the re- 
sults are expected to remain qualitatively valid for 
conventional optical fibers. After substituting Eqs. (1) 
and ( 2 ) into Maxwell’s equations, and making the slowly- 
varying-envelope and rotating-wave approximations, we 
obtain the following equations for the slowly varying com- 
plex amplitudes A and B (Ref. 4): 


8A 

i 1 I Bn fiA 


— 

= - B - - aA 5 - 


Sz 

2 2 2 dt 2 



+ (0+ iy)\A\ 2 A, 

(3) 

d B 

df 

= US - 1 )B - iAg, 

(4) 


frequency of the laser spectrum. The main assumptions 
in our model are the homogeneously broadened gain me- 
dium and the neglect of spontaneous emission. The 
former is not valid for all doped fibers, but for some types 
of glass hosts it is a reasonable assumption. 4 Since we 
are interested in deterministic instabilities, spontaneous 
emission can be neglected without loss of generality. 

There are two distinct origins of the nonlinear effects in 
Eqs. (3M5). The fiber nonlinearity y = n 2 w 0 /cA efr ac- 
counts for SPM effects induced by the host, where n 2 is 
the nonlinear refractive index (units m 2 /W), c is the speed 
of light in vacuum, and A efr is the effective fibercore area. 
For completeness we give the relation between n 2 and the 
nonlinear susceptibility x i3) (units of meters squared per 
volt squared) of the fiber 

3 

n 2 = “ Re [X' <3 Wo)]> <6) 

4 e 0 rc c 

where n is the background refractive index and e 0 is the 
permittivity of the vacuum. The dopant-induced nonlin- 
ear effects are governed by the saturation power P sat , de- 
fined as 

h 2 cn e(>A efr 

^sal “ ~ 7 ) 271 rp * ^ 

2^ TiT 2 

where h is Planck’s constant divided by 277 and g is the 
dipole moment of the atomic transition. Note that Eqs. 
(3)-(5) are written in the frame of reference moving with 
group velocity v g s which means that t - T 

- fiiz, where T is the time in the rest frame. By doing 
this, we eliminate the term p x idAfdT) from the left-hand 
side of Eq. (3). 

We now briefly discuss the relation of Eqs. (3M5) with 
the CGL model. 4,12 When the assumption is made that 
the population relaxation time T j is much longer than all 
other lifetimes, we can approximate the actual gain g by 
its steady-state value g 8 . This allows Eq. (4) to be ex- 
pressed in the Fourier domain as the well-known 
Lorentzian-shaped nonlinear susceptibility: 


d g , 

Tl dt = e ° “ g + Im(A B)/P (5) 

where g is the gain realized by pumping the dopants, a is 
the optical loss, /3 2 is the GVD coefficient of the host fiber, 
8 accounts for saturable host absorption, y is the fiber 
nonlinearity, 8 - (w 0 - w)T 2 is the scaled detuning be- 
tween the carrier frequency w 0 and the atomic resonance 
frequency w, g 0 is the unsaturated gain, and P Mt is the 
saturation power for the dopants modeled as a homoge- 
neously broadened two-level system. We have written 
Eqs. (3M5) in such a way that A has units of >/W, B has 
units of \fWL~ l , and g has units of L' 1 , where L is the 
length of either the amplifier or the laser cavity. 

It is important to note that Eqs. (3)-(5) are based on a 
traveling-wave description rather than a standing-wave 
approach that is employed in the conventional rate- 
equation analysis. Since we adopt a traveling-wave ap- 
proach, the optical field A (z, t ) in Eqs. (3)— (5) in principle 
can represent a very wide spectrum (or many longitudinal 
modes). The detuning parameter 8 is thus interpreted as 
the mismatch between the gain peak and the dominant 


A(Ao>) 1 - i(A(jjT 2 + 8) * 

where Aw = w - u> 0 is the detuning of the spectral com- 
ponent from the carrier frequency. In the CGL model, 
polarization equation ( 8 ) is approximated by a Taylor ex- 
pansion near Aw = 0 up to second order, leading to the 
parabolic-gain approximation that is reasonably accurate 
for small values of T 2 . In the time domain this corre- 
sponds to (generally complex) corrections Afi x and A(3 2 °f 
the inverse group velocity and the GVD coefficient / 3 2 
(Ref. 3): 

4W<I = ^ l9> 
-e.T, \ ‘ (10) 

(1 — id) 

From Eq. (9) we see that only at resonance (8 = 0 ) can 
the resulting pulse propagation equation be written in the 
reference frame moving with the new group velocity (/?j 
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+ A/?,)' 1 , because only at resonance is the correction 
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ppt 0 j *,v“ ^ s7,2 ’ thC resultin £ equations of the 
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( 12 ) 


t 0 het I n Sly ’ ( ! !) “ d (12) are ^ 00d approximations of 
the full model only if T l is long enough and T 2 is short 

enough. In this paper we explore the shortcomings of the 
CGL model for realistic fiber lasers and amplifiers and 
find interesting behavior outside the realm of the CGL 

Ztm r n0t t th , at the CGL model is onl y useful for 
amplifiers for which gain saturation can be neglected- 

otherwise, the gain dispersion b would be r dependent 
which seems impractical at best. 

The Maxwell-Bloch equations (3)-(5) can be applied to 
both amplifiers and unidirectional (e.g., ring) lasers. In 
the case of lasers, however, one should, in general solve a 
complicated boundary-value problem to account for the lo- 
calized losses at the cavity mirrors, a task that requires a 
numerical approach. In this paper we adopt the mean- 
intensity approximation by replacing the localized mirror 
osses with a distributed loss incorporated in the total op- 
tical loss a. In the case of amplifiers, such a mean- 
requirement is not valid: the intensity is 

strongiy z dependent. Because the steady states are so 
different for lasers and amplifiers, a modulation stability 

™ iis ■» <*■ 
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We consider an amplifier (or absorber) of length L with an 
p power P 0 at z - 0. We first find the time- 
independent (steady-state) solution of Eqs. (3)-(5). For- 
mally, it can be written as 

•A,(z) = [P^aJJ^exp^.U)], 


B,(z) = 


A $(z)g.(z) 


8 + i 


E>(z) - g 0 


1 + 


p a<z) 

P*( 1 + <S 2 ) 


(13) 

(14) 

(15) 




Using the imaginary part of Eq. (3), we can write the 
phase profile <p,( z ) ln terms of the power profile P A (z): 


<P„(z) 


= rf’dz' 

JO 


Pa( z ') + 
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1 + J 2 


/ dz'g,(z'). 
Jo 


(16) 


This differential equation can be solved, resulting in 
following transcendental equation for P A 

20P„,(1 + S 2 )z 


_/ j 1 1 " C\ P A (z) 

C_ '[ P 0 j 

+ - c J Pa<2) ~ c ^ p “ t(1 + s2) 


From the real part of Eq. (3), and using Eq. (15) one finds 

fte followtag JUtottU equation hr & 

profile f(z ) . P A (z)/[P ut ( 1 + s 2 )): P 
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where the coefficients C and C- are given by 
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In the absence of saturable absorption ( 0 = 0) the sol 

tion of Eq. (17) is implicitly given by 
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which, in the absence of optical loss (a = 0), reduces 


In 


PaU) 


P a(Z ) P 0 ggZ 


Bat 


1 + s 2 ' 


^ £ We T ofile Pa{Z ) ,S f0Und - the gain prof 
fl3)-(16) d Ule P ° lanzation P rofile B >(z) follow from Ec 

To study the onset of MI, we follow a standa: 
approach by considering the linear stability of tl 
steady-state (cw) solution given above. Considerir 
small perturbations u, v f p a and x frnm * 

defined as P ' 9 ’ * fr ° m the w sta ' 
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A(z, t) = [(P 0 r + u{z, t) + iviz, t)] 


Pa(z) 


P o 
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C23) 
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\Pa(z)} 


m 


exp[i<p a (z)], 
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g(z , t) = [g 0 + x(z, t)]\ 1 


+ 


Pa(z) 


P > aid + *) 


(25) 


and linearizing Eqs. (3)-(5) in u, v t p t q, and x, we solve 
the resulting five linear equations in Fourier space by in- 
troducing 


y(z, t) = y Q exp 


r 

- 


i dzK{z ) - iilt 


y = u y v, p, q y x y (26) 

where y 0 is the initial amplitude, fl is the frequency, and 
K(z) is the local wave number of the perturbation. The 
resulting dispersion relation for arbitrary detuning S is 
discussed for fiber lasers in Section 5. The dispersion re- 
lation at resonance {S = 0) is given by 


{[2 iK(z) + g s (z) - 40P A (z)][2iK(z) + g 9 (z)} 

+ /J 2 2 n 2 [o 2 + S gn(^ 2 )n c 2 (z)]}(i - ifiT 2 ) 

x [(i - inT’jKi - inr 2 ) + 7u)] 

~ g s (z)[2iK{z) + g,( 2 )](l - inr 2 ) 

x [1 - HIT, - 7(2)] -g,(z) 

x [2 iK(z) + g a (z ) - 4 0P A (z)] 

x [d - inr,)(i - inr 2 ) + /(*>] 

+ g,Hz )[ i - iar, - 7 ( 2 )] = o, (27) 

where I(z) = P A (z)/P Mt , 8gn(ft) = ±1, and n r ( 2 ) 
= [4yP i ,(2)/ 1>3 2 | J 172 is the critical frequency, i.e., the 
maximum frequency for which MI is found to occur in the 
case of anomalous dispersion in a passive fiber. 4 Before 
we examine the implications of Eq. (27) in various re- 
gimes of parameter space, we note that the Rabi fre- 
quency is somewhat hidden: 


^Rabi(^) s 


Pa(z) 


P MtP 1^2 


1/2 


(28) 


The imaginary part of K(z) determines the local gain 
experienced by the perturbation. It is useful to define 
the total integrated gain at frequency O as 11 


h(n) = 


-2 


f 


dz Im[tf(n, z)] y 


(29) 


where the factor 2 converts /t(H) to power gain. MI oc- 
curs whenever the wave number K has a negative imagi- 
nary part. In the case of an amplifier, this means the 


perturbation grows faster than the steady-state power, 
whereas for an absorber it means that the perturbation 
dampens less quickly. Dispersion equation (27) reduces 
to the previously reported ones in the appropriate limits. 
In the absence of saturable host absorption (0 = 0), the 
dispersion relation from Ref. 11 is obtained in the limit of 
large T x and short T 2 . 

We now consider the occurrence of MI in various re- 
gimes of parameter space for both amplifiers and absorb- 
ers. Because our model has no restrictions with respect 
to the magnitude of the lifetimes T x and T 2 , we can ex- 
plore MI in regimes where the CGL model has no validity. 
For simplicity, we only consider the local perturbation 
gain because the integration in Eq. (29) can be performed 
analytically in a few limiting cases only. 3 We also ignore 
the possibility of saturable host absorption since two- 
photon absorption is relatively weak in silica fibers, and 
other sources of saturable nonlinearity are rarely present 
in amplifiers. When we discuss MI in lasers, we show 
how even relatively small amounts of saturable host ab- 
sorption can affect the MI drastically. 

We start by investigating the effect of the magnitude of 
the dipole lifetime T 2 . For most fiber amplifiers, T 2 is 
estimated to be near 100 fs, corresponding to a wide gain 
spectrum. Because is usually in the range 0,1-10 ms, 
the CGL equation is expected to be a good approximation. 
However, by cooling the fiber, the polarization dephasing 
process can be slowed down substantially, making values 
of T 2 ~ 10 ps readily attainable. 13 

In Fig. 1 we show for various values of T 2 in the range 
0.1-10 ps the MI spectrum for a typical fiber amplifier 
with a 30-dB gain, i.e., exp(g 0 L) = 1000. All other pa- 
rameters are given in the caption. Note that the satura- 
tion power P aat is inversely proportional to T 2 [Eq. (7)]. 
When Eq. (7) is satisfied for each value of T 2 , the Rabi 
frequency remains a constant for all curves: 

= 1.29 X 10" 3 n c . When the dephasing time T 2 is in- 
creased, two trends are observed. 

First, as can be seen in Fig. 1, increasing T 2 leads to a 
shrinkage of the MI bandwidth, whereas the maximum 



• 0 0.5 1.0 US 


Frequency [O/Q, ] 

Fig. 1. Modulation instability spectrum for an erbium-doped fi- 
ber amplifier at various values of T 2 (indicated in the figure). 
Parameters are g Q - 6.91 L“\ P 0 = 1 mW, T Y = 0.1 ms, 0 2 
= -20 ps 2 /L, r=3r‘ L" 1 , and ? ul =lmW when T 2 
= 0.1 ps. For undoped fibers, MI occurs up to 2tt 
= 3.9 GHz 
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Fig. 2. MI spectrum for the amplifier of Fig. 1, for even longer 
dephasing times T 2 (indicated in the figure). In the range 21 
< T 2 < 80 ps, MI is totally quenched. When T 2 approaches 
100 ps, the MI spectrum starts to show a narrow, weak peak 
around the Rabi frequency = 0.0013. 

gain within the bandwidth decreases rapidly. When T 2 

- 100 fs ( P A /P sat = 0.01), MI occurs for frequencies up 
to n c /27T = 3.9 GHz, while the peak MI strength is found 
near ft = fl c / v2. The peak strength is very close to the 
analytical value 2 yP A , which is found in the CGL limit. 11 
Already when T 2 = 2 ps (P A /P aat — 0.2), the frequency 
band where MI occurs has shrunk -40 9c t and near T 2 

— 8 ps ( P A /P s at = 0.8), MI has almost ceased to occur at 
all. Near T 2 — 21 ps (P A /P 8at = 2.1), the MI band van- 
ishes completely. Long before that happens, the strength 
of MI is so weak that it is doubtful whether it can be ob- 
served in a single-pass amplifier. 

Second, during this MI spectrum shrinkage, another 
phenomenon is occurring that is directly caused by the 
two-level system since it involves frequencies close to 
^Rabu as is shown in Fig. 2. Near the Rabi frequency a 
secondary, weak maximum in MI strength begins to form 
for T 2 > 20 ps. This maximum becomes positive near 
T 2 = 80 ps (P A /P Bat = 8) and grows with T 2 . When T 2 
is increased further, the MI spectrum slowly returns to its 
original width and strength (out of scale in Fig. 2). Near 
T 2 = 11.5 ns (P A fP gat = 1150), the MI spectrum shows 
again positive MI gain around fl c , while maintaining a 
narrow (but weak) peak close to fl^bi . At the highly im- 
probable value of T 2 - 1 M s (P A /P ial = 10 5 ) the MI spec- 
trum is very close to the one at T 2 = 100 fs, and we have 
come full circle. 

Thus we find four regimes of T 2 \ in the first regime 
( 100 fs < T 2 < 21 ps), increasing T 2 leads to a total 
quenching of MI. In the second regime (21 ps < T 2 

< 80 ps). no MI occurs, but the gain around the Rabi fre- 
quency is growing. In the third regime (80ps< T 2 

< 11.5 ns), more and more MI occurs around the Rabi 
frequency, while the gain around fl r is growing toward a 
positive value again. In the fourth regime, approaching 
the long T 2 limit (1 M s < T 2 < =c), the MI spectrum re- 
covers fully to its original (small T 2 ) form. The bound- 
aries between these regimes are, of course, strongly de- 
pendent on the power level P A . For higher power levels 
these boundaries rapidly decrease. 


We further note that the MI band near 11,. is insensit 
to changes in T x , as long as it is accompanied by a chat 
in the saturation power P sat according to Eq. (7). Hi 
ever, if we keep the saturation power constant up 
changing (this can be done by adjusting the dipole n 
ment n), decreasing T , leads to a stabilization of t 
lower frequencies and eventually a reduction of MI al 
gether. 

We emphasize that the narrow MI peak around ll Rab 
so weak that it is questionable whether it can be observ 
in an amplifier. In the case of a laser, however, such 
weak gain may build to a substantial instability ov 
many round trips, as we discuss in the next sectu 
Since it is not common to use an amplifier in the high 
saturated regime, the emergence of the narrow MI bai 
near the Rabi frequency is not very practical. Note, ho\ 
ever, that this narrow MI band near the Rabi frequen. 
does not depend on the sign of (3 2 ; in both normal- ar 
anomalous-dispersion regimes, this instability emerges ; 
relatively high values of T 2 . 

Apart from this new (and for realistic systems, e> 
tremely weak) instability, the full Maxwell-Bloch modi 
agrees with the CGL model qualitatively rather well. C 
course, the quantitative differences become larger as th 
approximations leading to the CGL model (large T, an 
short T 2 ) become more and more inappropriate. In th 
next section we find that for lasers the situation can b 
very different. 


muuulahun instability in fiber 

LASERS AT RESONANCE 


Equations (3)-(5) also describe the optical field and tht 
gain in a laser, when one assumes that all losses can bt 
thought of as being distributed along the cavity. Then 
the steady-state solution is characterized by a 
2 -independent power P 0 and gain g, , and can be written 
as 


A s (z) - (P 0 ) 1/2 exp[t>,(z)], (30 1 

S - i 

B ’ iz) = J~ + J2 (31) 


8 s 8 o 


1 + 


P«*( 1 + s 2 ) 


(32) 


Again, from the real and the imaginary part of Eq. (3), the 
following expressions for the laser power P 0 and the 
phase profile <p s (z ) are obtained: 


g s - (a - 2<9P 0 )(1 + S 2 ), (33) 


d<P s 

IT = yP ° + 


i/ sg* 

2 \ 1 + S 2 



(34) 


Since Eq. (33) is quadratic in P 0 [with use of Eq. (32)], in 
principle, two values for the laser power are found. One 
of these is not physical and corresponds to the antilaser, 
which is characterized by a huge gain and almost zero 
power. 

Similar to the amplifier case, we consider small pertur- 
bations u, v, p, q , and x from the cw state, defined as 
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A(z, t) = [(P 0 ) m + u(z , t) 

+ iv(z , 0]exp[i<p,(z)], 


(35) 


B(z, t) = 


S — i 

1 + < 5 2 


[( Pq ) 1 ' 2 


+ P(z, 


t ) 


+ iq(z, OJgjexptiv^z)], 

P o 


g(z, t ) = [g 0 + *(*. <)] 


1 + 


P« td + *) 


(36) 


, (37) 


and linearizing Eqs. (3)-(5) in u, v, p, q , and x, we solve 
the resulting five linear equations in the Fourier space by 
introducing 


y(z, t) = y 0 exp [i(Kz - lit)], 

y = u, v, p, q, x, (38) 

where y 0 denotes the initial amplitude of the perturba- 
tion. Note that because because both laser power P 0 and 
g, are z independent, the wave number K is also z inde- 
pendent. At resonance (<$ = 0), the resulting dispersion 
relation reads 

{[2iA" + a - 60P 0 ](2iK + a - 2 OP 0 ) 

+ /? 2 2 n 2 [H 2 + sgn(/? 2 )ll c 2 ]}(l - inr 2 ) 
x [(i - inr,)(i - inr 2 ) + / 0 ] 

- gs (2iK + a - 28P 0 )(1 - ifir 2 ) 

x [1 — iHT 1 ! — / 0 ] - g,(2iK + a - 60Po) 

x [d - inrod - inr 2 ) + / 0 ] 

+ g s \ 1 - iClT x - 7 0 ) = 0. (39) 

Here, J 0 = sgn(ft) = ±1, and H c 

= (4yP 0 /|/3 2 |) 1/2 is the critical frequency, i.e., the maxi- 
mum frequency for which MI is found in the case of 
anomalous dispersion in a passive fiber. 4 This dispersion 
relation is identical to Eq. (27) when one replaces K , P 0 , 
andg, by their 2 -dependent counterparts, and Eq. (33) is 
used. 

Dispersion relation (39) reduces to the one previously 
reported by Chen et al., 12 who employ the CGL model, in 
the appropriate limit. 

Before we proceed with examining the implications of 
Eq. (39) in various regimes of parameter space, we note 
that the Rabi frequency is now given by 

/ P \ 1/2 

• (40) 

Merely comparing the relative strengths of the critical 
frequency fl c with the Rabi frequency H^bi does not pro- 
vide much information about the effect of atomic coher- 
ence on ML The interaction between the fiber nonlinear- 
ity, the GVD, and the two-level system is much more 
involved. 

In contrast with the amplifier case described in the pre- 
vious section, lasers generally operate in the heavily satu- 
rated regime. This means that the instability near the 
Rabi frequency is now more likely to play a significant 
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net Ml gain [-2 lm(K)] 


Fig. 3. MI analysis for a figure-eight laser. Solid curves indi- 
cate the results of the full model, while dashed curves show those 
of the CGL model. Top figure shows the net MI gain spectra, 
while the bottom figure shows the corresponding trajectory of the 
eigenvalue K on the complex plane. Parameters are a 
= 0.4 L" 1 , g 0 = 6 L"\ = -0.09 ps 2 6= 0.1 W 1 L~\ y 

= 0.008 W _1 L _1 , r 2 = 1.27 ps, 7\ = 10 8 ps, and P Mt = 10 
mW. 


role. Furthermore, many fiber laser systems, e.g., a 
figure-eight laser, contain an effective saturable absorber 
that causes mode locking. 

We first compare the predictions of Eq. (39) with the 
CGL-based expression. 12 To facilitate comparison, we 
used the same parameters as in Ref. 12. In Fig. 3 we 
show the differences for the case of a figure-eight laser. 
Although the trajectories of K in the complex plane as a 
function of frequency Cl are quite different for a figure- 
eight laser, the resulting net MI gain spectra agree quite 
well, at least in the central region. The frequency range 
over which positive net MI gain occurs is underestimated 
by 10% by the CGL model. Both models show vanishing 
gain at 100 kHz [indicated by the vertical line at ft — 0 in 
Fig. 3(a)], whereas the frequency with highest gain is 
near 200 kHz. The Rabi frequency in this case is 
= 55 MHz, and the_ critical frequency is n c — 37 GHz. 
They differ by almost three orders of magnitude. Even 
so, MI occurs for frequencies almost twice as large as H c . 

The results for the dye-laser parameters are shown in 
Fig. 4. Our model predicts that MI occurs in a narrow 
band near 30 GHz, whereas the CGL model predicts no 
instability at all! The Rabi frequency and the criti- 
cal frequency fl c are both close to 24 GHz, which explains 
why the interaction between the fiber and the two-level 
system is so highly nonlinear. 

Here we find the first meaningful qualitative difference 
between the full Maxwell -Bloch model and the CGL 
model. Not surprisingly, the laser power in Fig. 4 is —60 
times the saturation power P Mt , which makes the Rabi 
frequency of the same order as the critical frequency fl c . 
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So for the dye laser of Fig. 4 the interaction between the 
two-level system on the one side and the GVD and the 
SPM on the other side cannot be described within a 
parabolic-gain approximation. 

The pump value g 0 at which the cw state loses its sta- 
bility is often classified as the second threshold, 1 as it an- 
nounces the onset of unstable behavior. Similarly, we 
can identify the MI threshold as the gain above which MI 
occurs. At this threshold, MI occurs only at the fre- 
quency corresponding to the peak gain in the MI spec- 
trum, which can be compared with the frequency with 
which perturbations grow at a Hopf bifurcation. In Figs. 
5 and 6 we show the dependence of the MI threshold as a 
function of 6 for the case of the figure-eight laser (Fig. 5) 
and the dye laser (Fig. 6). The effect of saturable absorp- 
tion is very dramatic in the case of the figure-eight laser 


i 


<u 

c 



nJZn [THzJ 



Fig. 4. MI analysis for a dye laser. Similar as in Fig 3 except 
for the parameters: a = 0.1 L~\ g 0 = 3 L' 1 , /3 2 = -0.09 ps 2 

L-, 6 = 0.001 W*> L", y = 0.008 W' 1 L' 1 , T 2 = 2.45 ps T, 
= 10 3 ps, and P ul = 1 m W. 



^ threshold as a function of saturable absorption 6 for 
the fiber laser of Fig. 3. 



Fig. 6. MI threshold as a function of saturable absorption 6 for 
the dye laser of Fig. 4. 



Fig. 7. Comparison of net MI gam spectra in the absence of 
saturable absorption (0=0), for the cases of normal and anoma- 
lous dispersion. Other parameters of the fiber laser are the 
same as in Fig. 3. 


parameters: when <9is larger than 10“ 6 W' 1 L“ 1 . MI oc- 
curs immediately after the first (lasing) threshold. A 
qualitatively similar dependence is found for the dye la- 
ser, where the MI threshold gain decreases from ~ 9 to 
1.5 at 6 = 0.1. This feature explains why a relatively 
weak saturable absorber can lead to passive mode lock- 
ing. Although the presence of a saturable absorber is 
evidently very useful for the generation of mode-locked 
pulses, it somewhat obscures our investigation of the in- 
teraction between the two-level system and the fiber non- 
linearity and dispersion. This explains why for both the 
figure-eight laser and the dye laser of Figs. 3 and 4, the 
results hardly change if we consider normal dispersion. 

We therefore examine the case 0 = 0, so that we can 
consider the sole interaction between the fiber nonlinear- 
ity and the two-level system occurring in the absence of 
saturable absorption. The interesting question is 
whether atomic coherence can lead to MI in the normal- 
dispersion regime of the fiber. According to the CGL 
model, this is not possible. 

In Fig. 7 we show, again for a fiber laser, but without 
saturable absorption (ring-cavity instead of figure-eight 
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geometry), the net MI gain spectra for p 2 
= ±0.09 ps 2 L _1 . For anomalous dispersion, we find an 
approximately 4-GHz-wide MI band centered near 2.5 
GHz and a much narrower and much weaker MI band 
centered near 50 MHz. While the MI band near 2.5 GHz 
vanishes in the case of normal dispersion, the narrow low- 
frequency band survives. So, contrary to what the CGL 
model predicts, the presence of dopants can cause MI in- 
stability in the normal-dispersion regime of the fiber. 
Furthermore, in contrast to the amplifier case (Section 2), 
any positive value of —2 Im(if ) should be taken seriously, 
since, in a laser, even the smallest growth of a perturba- 
tion may cause a significant change in the output signal 
after many round trips in the cavity. 

Figure 8 shows the peculiar dependence of this new MI 
at normal dispersion when the population relaxation time 
Ti is decreased from 10 /is to 1.375 /is. Upon decreasing 
T x , the new MI band initially grows stronger, while shift- 
ing to higher frequencies. Decreasing T Y further causes 
the band to weaken and finally to vanish abruptly at 
~1.35 /is. 

We stress that our results indicate that a fiber ring la- 
ser, operating in the normal-dispersion regime, may show 
unstable behavior at high pump levels, even in absence of 
additional saturable absorbing mechanisms. At reso- 
nance, the strength of the MI in the normal-dispersion re- 
gime is rather weak, which would imply that the instabil- 
ity needs to build up during many round trips in the ring 
laser. In the next section we discuss the effect of detun- 
ing on the strength and the nature of this new instability. 


5. MODULATION INSTABILITY IN 
DETUNED FIBER LASERS 

With various experimental techniques, e.g., through the 
use of gratings, one can force a fiber laser to operate away 
from the gain peak. In our theory this means that we 
have to deal with the effect of detuning 6. Recall that 
there is no CGL version for the detuned case, as the group 
velocity becomes complex in the parabolic-gain approxi- 
mation. For arbitrary detuning, dispersion relation Eq. 
(39) reads 
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Fig. 8. New MI band at normal dispersion as a function of popu- 
lation relaxation time T x (indicated). Other parameters of the 
fiber laser are the same as those in Fig. 3. 


where 7 0 = Equation (39) is recovered by put- 

ting S = 0 and g s = a - 2 0P O . Equation (41) can be 
applied as well for amplifiers by treating P 0 , g s , and K as 
z -dependent quantities. 

We now use Eq. j41) to investigate the effect of detun- 
ing on the MI spectra shown in Fig. 7. The introduction 
of <5 into the problem makes the situation even more com- 
plex. Instead of only two frequencies, i.e., the critical fre- 
quency fl c and the Rabi frequency the problem 

now is governed by the interaction of three frequencies. 
In Figs. 9 and 10 we show the effect of detuning on the 
bandwidth and the strength of MI for normal and anoma- 
lous dispersion, respectively. Clearly, small detunings 
have a large effect on the occurrence of MI, and the sign of 
the detuning also matters. This spectral asymmetry is 
due to the fiber host nonlinearities. When GVD and 
SPM are absent, Eq. (41) is symmetric in detuning S. In 
the anomalous-dispersion regime (Fig. 10), a small value 
of the detuning connects the two MI bands, one owing to 
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Fig. 9. Effect of detuning on the new MI in the normal- 
dispersion regime. Parameters identical to those in Fig. 7. 



Fig. 10. Similar to Fig. 9, except that the laser now operates in 
the anomalous-dispersion regime. 

the passive fiber MI and the other owing to two-level dy- 
namics. Upon increasing the absolute value of the de- 
tuning, the MI bandwidth and the strength increase loga- 
rithmically. At large detuning (<5 = 0.1), there is no 
distinction between the normal- and anomalous- 
dispersion case. For large detuning the instability is ap- 
parently dominated by the two-level dynamics. 


6. CONCLUSIONS 

We have analyzed the occurrence of modulation instabil- 
ity (MI) in fiber lasers and amplifiers by considering the 
self-phase modulation, group-velocity dispersion, and the 
saturable host absorption. The gain spectrum has been 
fully considered, in contrast to the parabolic-gain approxi- 
mation employed in the complex Ginzburg-Landau 
(CGL) model. We have derived analytical expressions for 
the MI dispersion Kdl) that naturally reduce to previ- 
ously reported research for both lasers and amplifiers. 

For amplifiers, operating not too heavily saturated and 
in absence of saturable absorption, no qualitative differ- 
ences with the CGL description are found, even in re- 
gimes where the basic approximations of that model are 
violated. Quantitatively, however, the differences can be 
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quite substantial. We show that by cooling the fiber am- 
plifier and thereby increasing the dipole dephasing time, 
the occurrence of MI can be quenched. For heavily satu- 
rated amplifiers, we find a new instability located in a 
narrow frequency band around the Rabi frequency. The 
CGL model does not predict such an instability. The 
strength of this new instability is very small, and it is 
questionable that its effect can be detected in a single- 
pass amplifier. 

In lasers, a different picture emerges, since any grow- 
ing perturbation may build up over many round trips 
within the laser cavity. Furthermore, the presence of a 
weakly saturable absorbing mechanism is shown to 
greatly enhance the instability. We compare our results 
with those of Chen et a/., 12 who used the CGL model to 
investigate MI in a dye laser and a figure-eight laser. 12 
Our results for the figure-eight laser agree rather well, 
whereas we find disagreement for the dye laser, which in 
our model is predicted to have an instability of —30 GHz. 
Further indication that the CGL model should be used 
with caution is given when systems without saturable ab- 
sorption are studied: for a fiber ring laser operating in 

the normal-dispersion regime, a narrow MI band of low 
(—50- MHz) frequencies is found, which is not predicted 
by the CGL model. This may explain the self-starting of 
mode-locked Nd-doped fiber lasers. 

The effect of detuning on the strength and the band- 
width of the new instability can be substantial, since non- 
zero detuning effectively introduces a new frequency into 
the problem. Even for a relatively small detuning, the 
strength and the bandwidth of MI increase logarithmi- 
cally, whereas the difference between normal and anoma- 
lous dispersion becomes smaller. The fiber nonlineari- 
ties cause the MI spectrum to become asymmetric with 
respect to detuning. 
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Absolute Instabilities in Lasers with 
Host-Induced Nonlinearities and Dispersion 

Guido H. M. van Tartwijk and Govind P. Agrawal, Fellow, IEEE 


Abstract — We analyze the occurrence of absolute instabilities 
in lasers that contain a dispersive host material with third- 
order nonlinearities. Starting from the Maxwell-Bloch equations, 
we derive general multimode equations to distinguish between 
convective and absolute instabilities. We find that both self-phase 
modulation and intensity -dependent absorption can dramatically 
affect the absolute stability of such lasers. In particular, the self- 
pulsing threshold (the so-called second laser threshold) can occur 
at lew times the first laser threshold even in good-cavity lasers for 
which no self-pulsing occurs in the absence of intensity -dependent 
absorption. 

Index Terms — Laser stability, nonlinear optics, optical fiber 
lasers, optical Kerr effect, optical pulse generation, optical prop- 
agation in dispersive media. 

I. Introduction 

A LMOST immediately after the advent of the laser, it 
was recognized that laser output C3n become unstable, 
resulting in irregular power spikes even at a constant pumping 
level |lj. Over the last 30 years or so. laser instabilities 
have been studied extensively both from the fundamental and 
applied viewpoints [2], [3]. The fundamental studies have led 
to the flourishing field of optical chaos. On the applied side, 
the development of techniques for controlling chaos are being 
used to make lasers tailored for specific applications (high 
power, short pulses, clean far field, etc.). 

Since deterministic chaos is studied in a wide variety 7 of 
disciplines, the understanding of laser instabilities can be 
improved by referring to plasma and fluid instabilities that 
have been studied for a long time. A famous example is 
prov ided by the Lorenz-Haken equations which are named 
after the fluid dynamicist Lorenz and the laser theorist Haken 
[3]. [4]. In fluid dynamics, instabilities are categorized into 
two types: convective and absolute [5]. Convective instabilities 
are characterized by the growth of localized perturbations 
upon propagation inside a nonlinear medium, while absolute 
instabilities exhibit purely temporal dynamics. Absolute laser 
instabilities have been studied for more than 30 years. The 
Lorenz- Haken equations describe the dynamics of a homo- 
geneously broadened gain medium in a unidirectional ring- 
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cavity. Although rarely stated explicitly, the Lorenz-Haken 
equations can only show absolute instabilities. The fundamen- 
tal concepts such as second laser threshold self-pulsing, Hopf 
bifurcation, and different routes to chaos are all formulated 
within the context of absolute laser instabilities [3]. 

In the last 15 years or so. new laser systems have been 
designed that arc not easily modeled by the Lorenz-Haken 
equations. Examples of such lasers are fiber lasers and solid- 
state (e.g.. Ti : sapphire) lasers, which are capable of producing 
ultrashort optical pulses through passive mode locking while 
operating at a constant pump power. What these lasers have in 
common is that the gain is provided by atoms or ions doped 
inside a host material. As a result, the cavity' contains not 
only a gain element but also other nonlinear elements, which 
are responsible for nonlinear processes such as self-phase 
modulation (SPM) and intensity- dependent absorption (EDA) 

[6] , Also, group-velocity dispersion (GVD) of the host medium 
plays a nonnegligible role. Because of the dispersive and 
nonlinear effects, evolution of the optical field over a single 
round trip must be considered, contrary to the Lorenz-Haken 
model in which such effects are ignored. This means that the 
convective nature of any instability must be considered while 
discussing instabilities for such lasers. 

A well-known example of a convective instability occurs 
in nonlinear fiber optics [6]. Optical fibers, without any 
gain dement and without any longitudinal resonances (no 
cavity), show a convective instability known as the modulation 
instability. When the power of a CW optical beam becomes 
sufficiently large, the combination of SPM and anomalous 
GVD causes the CW beam to break up spontaneously into 
a pulse train (and eventually into optical solitons) whose rep- 
etition rate depends on the fiber parameters. Mathematically, 
a linear stability analysis shows that perturbations of the form 
exp[—?(D/ — K z )\ grow exponentially as exp(pz) with a 
growth rale g — — Im(/<') that depends on the frequency 
of perturbation. The repetition rate of the resulting pulses 
corresponds to the frequency D for which the growth rate g 
is maximum. 

Adding gain to the system, e.g., by doping the fiber with 
rare-earth ions and pumping it optically, can affect consider- 
ably the conditions under which modulation instability arises 

[7] . The instability, however, remains convective in nature. 
When such a host material (with or without gain) is put into a 
cavity, the resulting boundary conditions at the cavity mirrors 
can change the nature of the instability from convective to 
absolute. Feedback is a necessary ingredient for absolute 
instabilities to occur. A well-known example is the Dceda 
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instability [8], which arises when a Kerr medium is placed in 
a unidirectional ring cavity. Even without gain and dispersion, 
the feedback mechanism provided by the cavity results in an 
absolute instability. 

In this paper, we discuss under what conditions a convective 
instability becomes absolute in a laser. In Section IL we 
derive, starting from the Maxwell-Bloch equations, a set of 
Lorenz -Haken- type multimode equations capable of describ- 
ing the temporal evolution of a laser whose cavity includes 
optical elements exhibiting dispersion and nonlinearities. The 
usefulness of this new set of equations is illustrated in Section 
III by considering a relatively simple case of a single-mode 
laser. We discuss the stability of that mode as a function of 
host nonlinearities. 


II. Maxwell-Bloch Equations 

For definiteness, we focus on a fiber laser although the 
analysis can be applied to any solid-state laser with some 
modifications. Our starting point is a set of Maxwell-Bloch 
equations describing the propagation of optical fields in an 
optica] fiber, doped with rare-earth ions. We write the optical 
field £ and the dopant-induced polarization V as 

y , 2, t ) - ^xF{x, y)A(z, t ) exp[i(Ao2 - u>o<)] + c.c. 

" ( 1 ) 
V{r, y, 2, t) - y)B(z, t ) exp[i(tf 0 2 -wo<)] + c - c - 

( 2 ) 


where x is the polarization unit vector of light assumed 
to be linearly polarized along the x axis, F(x , y) is the 
transverse profile of the fundamental fiber mode, and Kq is 
the wavenumber corresponding to the optical frequency u?o. 
We assume that the field-polarization direction is preserved 
upon propagation. After substituting (1) and (2) in Maxwell’s 
equations, modeling dopants as a homogeneously broadened 
two-level system, and making use of the slowly varying 
envelope and rota ting-wave approximations, we obtain the fol- 
lowing equations for the slowly varying complex amplitudes 
A and B [6]: 


dA 

dz 


1 dA i 1 i() 2 d 2 A . 1a . 2a _ 

v g dt ~ 2 B 2 aA 2 dt 2 + A (3) 


dB 

— = — ( 1 - i6)B - iAg 
T x l ^=g c -g + \m(A*B)/P. 


(4) 

(5) 


where tj is the gain realized by pumping the dopants, a is the 
optical loss of the host fiber, 7\ is the population lifetime 
of the dopants. T 2 is the dipole-dephasing time, v g is the 
group velocity, is the GVT) coefficient of the host fiber, 
the complex parameter 7 accounts for the host nonlinearities 
responsible for SPM and IDA, b — (u>o — ua)T 2 is the scaled 
detuning between the optical frequency u>o and the atomic 
resonance frequency' <70 is the unsaturated gain, and P t is 
the saturation power for the dopants. We have written ( 3) — f5) 
in such a way that A has units of >/W, B has units of 
xAv • m l , and g has units of m l . 


The main assumptions in our model are the use of a 
homogeneously broadened gain medium and the neglect of 
the stochastic nature of spontaneous emission. The former is 
not valid for all dopants but is a reasonable assumption for 
many types of dopants [6]. The latter can be justified if one is 
interested only in deterministic instabilities. 

There are two distinct origins of the nonlinear effects in 
(3)-(5). The host nonlinearity 7 = 7' + 17" accounts for SPM 
and EQA effects induced by the silica fiber. The SPM effects 
are governed by y* = njwo/c^ff, where n 2 is the nonlinear- 
index coefficient, c is the speed of light in vacuum, and A^r is 
the effective mode area [6]. The effects of EDA arc accounted 
for by 7". When y N > 0, the loss in the cavity increases with 
intensity, modeling processes such as two-photon absorption 
[6]. In contrast, negative values for y /f imply a decrease in 
cavity losses with increasing intensity and model fast saturable 
absorption. The dopant-induced nonlinear effects are governed 
by the saturation power P t = h 2 cntoA^ff / (2fi 2 TiT 2 ), where 
h is Planck’s constant divided by 2?r, fi is the dipole moment 
of the atomic transition, and n is the background refractive 
index. 

The Maxwell-Bloch equations, together with the boundary 
conditions imposed by the laser cavity, provide the most 
general framework for studying laser instabilities. They are 
capable of handling both convective and absolute instabilities 
and can show transitions between them. However, their solu- 
tions require a numerical approach. Without host nonlinearitics 
(7 = 0) and without GVD (#> = 0), the steady-state solutions 
can be obtained, and their linear stability properties have been 
studied [9]. However, such an approach is quite cumbersome, 
and it is not easy to carry out the analysis after the inclusion 
of host nonlinearities and GVD. If one is interested only in 
absolute instabilities, an analytic approach can be developed, 
as discussed in the next section. 

in. Multimode Laser Equations 

Rather than solving (3)-(5) numerically, we make use of 
the fact that any cavity supports a set of longitudinal modes 
whose field distribution / m (z, <) reproduces itself after each 
round trip inside the cavity. These modes can be obtained by 
solving (3) with B — 0 (no gain in the fiber cavity) and using 
the appropriate boundary renditions at the cavity mirrors. For 
a high-Q laser cavity, one can distribute the mirror losses 
throughout the cavity and replace the fiber loss a in (3) with 
a T = a + ax/, where ax/ is the distributed mirror loss. The 
boundary condition then simply becomes A(L, t) = A(0, <), 
where L is the cavity length. For a Fabry-Perot cavity with 
mirror reflectivities Ri and R 2 , a x/ is given by 



The approximation that the localized mirror loss can be 
replaced by a distributed loss only holds for a high-Q laser 
cavity [2]. 

When dealing with a unidirectional ring laser without host 
dispersion and nonlinearitics, the form of f m (z, t) becomes 
simply exp[— iu/ m (< — z/v g )], where are the mode fre- 
quencies and the loss term has been ignored. In the presence of 


